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A hybrid transport approach for the bulk evolution of viscous QCD matter produced in ultra- 
relativistic heavy- ion collisions is presented. The expansion of the dense deconfined phase of the 
reaction is modeled with viscous hydrodynamics while the dilute late hadron gas stage is described 
microscopically by the Boltzmann equation. The advantages of such a hybrid approach lie in the 
improved capability of handling large dissipative corrections in the late dilute phase of the reaction, 
including a realistic treatment of the non-equilibrium hadronic chemistry and kinetic freeze-out. By 
varying the switching temperature at which the hydrodynamic output is converted to particles for 
further propagation with the Boltzmann cascade we test the ability of the macroscopic hydrody- 
namic approach to emulate the microscopic evolution during the hadronic stage and extract the 
temperature dependence of the effective shear viscosity of the hadron resonance gas produced in the 
collision. We find that the extracted values depend on the prior hydrodynamic history and hence 
do not represent fundamental transport properties of the hadron resonance gas. We conclude that 
viscous fluid dynamics does not provide a faithful description of hadron resonance gas dynamics with 
predictive power, and that both components of the hybrid approach are needed for a quantitative 
description of the flreball expansion and its freeze-out. 

PACS numbers: 25.75.-q, 12.38.Mh, 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 

The discovery that ultrarelativistic heavy-ion colhsions 
at the Relativistic Heavy-Ion Collider (RHIC) produce 
dense, color deconfined matter that thermalizes quickly 
into a quark-gluon plasma (QGP) [1] and subsequently 
evolves like an almost perfect liquid [2-5] with near- 
minimal viscosity [6, 7] has generated intense interest in 
a quantitative determination of the QGP transport prop- 
erties [8-12]. Since the QGP liquid stage is sandvifiched 
between an early pre-equilibrium and a final hadronic 
rescattering and decoupling stage, both of which have 
different transport properties, such a determination re- 
quires a complete dynamical description of all stages of 
the fireball expansion [11]. Stated differently, when one 
uses experimental observables that are sensitive to the 
transport properties of the expanding medium (for ex- 
ample, the elliptic flow V2 has been shown to be par- 
ticularly sensitive to shear viscosity [8-13]), contribu- 
tions from the pre-equilibrium and hadronic freeze-out 
stages to these observables must be accurately known 
for a precise determination of the QGP transport coef- 
ficients. Purely hydrodynamic calculations, that substi- 
tute hydrodynamic initial conditions for a full dynamical 
solution of the pre-equilibrium stage and replace the ki- 
netic hadron freeze-out process by a sudden transition 
from thermalized fluid to non-interacting particles using 
the Cooper-Frye prescription [14], require additional pa- 
rameters whose values may influence the extracted ther- 
mal and transport properties, and they will always be 
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plagued by irreducible systematic uncertainties in the ex- 
tracted values that result from the crude modeling of the 
early and late non-equilibrium stages. 

This insight is not new, and it has spurred the de- 
velopment of hybrid algorithms that describe different 
stages of the expansion with different tools for the last ten 
years. The main advantage of a hydrodynamic descrip- 
tion (where valid) over kinetic theory is its relative sim- 
plicity: only a few macroscopic fields (energy and baryon 
density, pressure and flow velocity) must be evolved in 
space-time, whereas the microscopic approach requires 
to follow the evolution of both the momenta and po- 
sitions of all particles. The phenomenological success 
of hydrodynamic modeling for heavy-ion collisions thus 
leads naturally to the concept of using a fluid dynam- 
ical description as the backbone of the complete evo- 
lution model, interfaced with computationally more de- 
manding microscopic algorithms to describe the early and 
late non-equilibrium stages. For the hadronic rescatter- 
ing phase several microcopic algorithms that solve cou- 
pled Boltzmann equations for the hadron distribution 
functions were developed and extensively tested in the 
1980s and '90s [15-19]. Hybrid codes that coupled an 
ideal fluid dynamical description of an expanding QGP 
to such hadronic rescattering codes and compared the re- 
sults with purely fluid dynamical calculations began to 
appear about ten years ago [20-23]. Here we present 
the first hybrid model that interfaces viscous relativistic 
hydrodynamics (specifically the (2-t-l)-dimensional algo- 
rithm VISH2+1 [10]) with a hadronic Boltzmann cascade 
(specifically the Monte Carlo algorithm UrQMD [18]), via 
the "hydro-to-OSCAR" converter H2O that converts hy- 
drodynamic output into hadrons, with positions and mo- 
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menta given in OSCAR format^ that can be read by UrQMD. 
We call this hybrid code VISHNU.^ In the present version 
of VISHNU the pre-equilibrium stage is not yet described 
dynamically, but continues to be replaced by initial con- 
ditions for the fluid dynamic evolution. 

An essential ingredient in VISHNU is the use of a state- 
of-the-art equation of state (EOS) for hot QCD matter, 
s95p-PCE [24, 25], which incorporates our best knowl- 
edge of the relation between pressure, energy and en- 
tropy density in the deconfined QGP stage from Lattice 
QCD [26] (see also [27]) and matches it to a realistic 
hadron resonance gas (HRG) at low temperatures, tak- 
ing into account that the abundances of stable hadrons 
(after strong decays of unstable resonances) are exper- 
imentally known to freeze out at Tchom ~ 165 MeV [28]. 
This requires the introduction of properly adjusted, tem- 
perature dependent non-equilibrium chemical potentials 
for all hadronic species below Tchcm [31-34]. Chemical 
freeze-out is immediate"^ and automatic in UrQMD, due to 
the rapid three-dimensional expansion of the fireball in 
the hadronic stage [29]. This implies that the hydrody- 
namic output that is fed into UrQMD must have the correct 
chemical composition, since otherwise the final hadron 
abundances disagree with experiment. If we decrease the 
switching temperature Tgw for the transition from hydro- 
dynamics to hadron cascade below the chemical freeze- 
out temperature Tchem, the hydrodynamic output thus 
has to reflect the correct partial chemical equilibrium 
(PCE) abundances that would, after resonance decays, 
produce the correct final yield ratios. This is also impor- 
tant for the elliptic flow: Although the non-equilibrium 
chemical potentials have a very small effect on the equa- 
tion of state p(e) [32], the distribution of the total mo- 
mentum anisotropy among hadron species strongly de- 
pends on the chemical composition of the hadronic sys- 
tem [32-35]. A comparison between chemical equilib- 
rium and partial chemical equilibrium EOS in ideal hy- 
drodynamics shows that the breaking of chemical equi- 
librium in the HRG can increase the differential elliptic 
flow V2{pt) for pions by 25% [32, 34]. 

VISHNU was developed to remove uncertainties from 
non-equilibrium dynamics during the late hadronic stage 
when extracting QGP transport coefficients from experi- 
mentally measured final hadron spectra. The main goals 
of the present article, however, are purely conceptual: 
by comparing the hadron spectra and differential elliptic 
flow from VISHNU with those from purely hydrodynamic 
simulations with different (temperature dependent) val- 
ues for the shear viscosity to entropy density ratio rj/s 
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^ "Viscous Israel-Stewart Hydrodynamics aNd UrQMD". 
^ It is a consequence of the smallness of particle-changing inelas- 
tic cross sections when compared with the much larger elastic 
and (resonance dominated) quasielastic scattering cross sections 
that keep the hadron gas close to thermal (but not chemical) 
equilibrium for a range of temperatures below Tchem [30] . 



during the HRG stage, we want to establish to what ex- 
tent the microscopic UrQMD dynamics of the HRG phase 
can be mimicked by an effective macroscopic calculation 
based on viscous fluid dynamics. By varying the switch- 
ing temperature between VISH2+1 and UrQMD wc explore 
the existence of a "switching window" , i.e. of a range of 
temperatures within which both VISH2+1 (with an appro- 
priate choice of transport coefficients) and UrQMD provide 
a valid description of the fireball evolution. Our goal here 
is not to compare VISHNU with experimental data from 
RHIC; this is done elsewhere [12]. Consequently, we ex- 
plore the sensitivity of our results to variations in the 
hydrodynamic initial conditions only to the extent that 
they affect the answer to the above questions. 

The paper is organized as follows. In Sec. II we briefly 
present the three components of VISHNU: the viscous hy- 
drodynamic algorithm VISH2+1, the hydro-to-micro con- 
verter H2O, and the hadron cascade UrQMD. In Sec. HI 
we compare basic observables (spectra and elliptic flow) 
from VISHNU and pure viscous hydrodynamics, using 
both chemical equilibrium and PCE equations of state 
as input. In Sec. IV we test the switching temperature 
dependence of the flnal spectra and elliptic flow obtained 
with VISHNU, extracting an effective viscosity {rj/s){T) 
for UrQMD under dynamical conditions provided by RHIC 
collisions. We interpret our flndings and summarize our 
results in Sec. V. Appendix A describes tests of the 
hydro-to- micro converter H2O and shows some results that 
demonstrate its accuracy. 

II. VISHNU: COUPLING VISCOUS 
HYDRODYNAMICS TO A HADRON CASCADE 

In this section we describe the structure of VISHNU, a 
hybrid code that couples the viscous hydrodynamic ex- 
pansion of the QGP stage in heavy-ion collisions to a 
microscopic kinetic evolution of the dilute late hadronic 
stage using a Boltzmann Monte Carlo approach. The 
hydrodynamic component allows for ideal fluid evolution 
in the limit of zero transport coefficients. We here in- 
clude only shear viscosity, neglecting bulk viscous contri- 
butions from the QGP and hadronization stages whose 
effects on spectra and elliptic flow are expected to be 
small [36, 37]."* We assume zero net baryon density ev- 
erywhere and thus do not follow explicitly the evolution 
of the baryon current. The latter can be easily included 
in future versions of the code. 

The viscous hydrodynamic and microscopic kinetic al- 
gorithms are interfaced with each other through the 
Monte Carlo event generator H2O that converts hydro- 
dynamic output on a hypersurface of constant temper- 
ature Tgw to particles, by sampling the corresponding 
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Cooper-Frye [14] phase-space distributions, including vis- 
cous corrections. This procedure requires switching tem- 
peratures Tsv, that are low enough for UrQMD to pro- 
vide a reliable description of the subsequent hadronic 
rescattering dynamics. The highest possible switching 
temperature is therefore Tgw = Tc where Tc denotes the 
(pseudo)critical temperature for the quark-hadron phase 
transition. Higher Tgw values would require a microscopic 
description of the hadronization process itself, including 
accompanying changes in vacuum structure, and knowl- 
edge of the effective particle degrees of freedom during 
this process. This is at present an unsolved problem. 

For comparison we also perform purely hydrodynamic 
simulations without the hadronic cascade, by running 
VISH2+1 to lower temperatures and decoupling directly 
into free-streaming particles, by using the Cooper-Frye 
prescription at Tdcc- These comparison runs are done 
with constant or temperature-dependent specific shear 
viscosities 11 /s in the hadronic phase. For the QGP phase 
we use constant rj/s, motivated the weak (logarithmic) 
temperature dependence of 77/s predicted by perturba- 
tive QCD in the weak-coupling limit [39] and the tem- 
perature independence of i]/s predicted by the AdS/CFT 
correspondence in the strong coupling limit [7].^ 

In the following subsections we discuss the components 
of VISHNU in more detail. 



from the optical Glauber model whose peak energy den- 
sity in central (6 = 0) Au-|-Au collisions is normalized 
to eo = e(r=0, tq; 6=0) = 30 GeV/fm^, with tq = 0.6 fm/c. 
This gives roughly the correct final charged hadron mul- 
tiplicity dNch/dy in central collisions, but does not cor- 
rectly reproduce its measured nonlinear dependence on 
the total number of wounded nucleons A^part in noncen- 
tral collisions. For the purposes of the present conceptual 
study this is not essential. 

In Sec. IV, in order to test the "universality" of the 
effective temperature-dependent (77/5) (T) for UrQMD ex- 
tracted by comparing VISHNU with VISH2+1, we also use 
Color Glass Condensate (CGC) motivated initial condi- 
tions, obtained by averaging a large number of fluctu- 
ating initial entropy density profiles computed with the 
f KLN Monte Carlo code® [40] after recentering them and 
aligning their major axes.^ The resulting smooth average 
entropy density is converted to an initial energy density 
profile using the EOS. This procedure accounts, in an 
average way, for event-by-cvent fluctuations in the initial 
source eccentricity, giving (depending on impact param- 
eter) 30 — 100% larger initial eccentricities than the 
optical Glauber model (see [41] for details). For both 
optical Glauber and fluctuating f KLN initial profiles we 
assume zero transverse flow at the beginning of the hy- 
drodynamic evolution at tq. 



A. Viscous hydrodynamics (VISH2+1) 

For the hydrodynamic stage we use VISH2+1, a (2+1)- 
dimensional viscous hydrodynamic code with longitudi- 
nal boost invariance [10]. The specific implementation 
used in the present work, including the equation of state 
s95p-PCE, is described in Sections II and III of Ref. [25] 
to which we refer the reader for technical details. For 
later reference we note that the relaxation time for the 
shear pressure tensor is set to = Zr] / {sT) . In some of 
the comparison runs we also employ the equation of state 
SM-EOS Q described in [10] which assumes a first-order 
quark-hadron phase transition and chemical equilibrium 
(CE) in the hadronic phase. The comparison between 
pure hydrodynamic simulations with s95p-PCE and SM- 
EOS Q (CE) emphasizes the effects arising from the 
breaking of chemical equilibrium in the hadronic phase 
at temperatures below Tchcm = 165 MeV. 

Different from Ref. [25], we here use as default ini- 
tial conditions an initial energy density profile taken to 
be 100% proportional to the wounded nucleon density 



^ We acknowledge that near Tc QCD is not a conformal field theory 
and that therefore in the temperature region accessible at RHIC 
(Tc <T<2Tc) r)/s might very well feature a stronger tempera- 
ture dependence than predicted by both perturbative QCD and 
the strong coupling limit for conformal field theories. Heavy- 
ion experiments at the LHC are expected to shed light on this 
question. 



B. Hydro-to-micro converter (H2O) 

To convert the hydrodynamic output into particles 
that can then be further propagated with the hadron 
cascade UrQMD we first use the AZHYDRO algorithm to find 
an isothermal freeze-out surface S(a;) of temperature Tgw 
and then calculate the hadron spectra with the Cooper- 
Frye formula [14] (see [10] for details). H2O is a Monte 
Carlo event generator that samples tiles on E and gen- 
erates the position x and momentum p of a particle of 
species i with a probability derived from the differential 
Cooper-Frye formula: 

= j£:^p-d'a{x)Mx,p) (1) 
= J^^P ■ d?(^{.x) [/oq,j(a;,p) + 5f^{x,p)] . 

Here gi is the degeneracy factor for hadron species i, 
and (Pafj,{x) is the outward-pointing surface normal vec- 
tor of the selected tile at point x on the surface E. 
fi{x,p) = fcq,i{x,p) + 5fi{x,p) is the local distribution 
function for hadron species i, consisting of a local cqui- 
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librium part (here with /i = 0) 

fcq.i{x,p) = gp.„(^)/T(a;) _^ ^ 

and a (small) deviation Sfi from local equilibrium due 
to shear viscous effects. For 6f we make the quadratic 
ansatz (see [8, 42] for a discussion of this and other pos- 
sibilities) 

5f{.,p)=Up, X) (lT/cq(p, X)) 2r2|l7(r(xHt(x)) 

(the upper (lower) sign is for fermions (bosons)). 5f is 
proportional to the shear viscous pressure tensor tt'^''{x) 
on the freeze-out surface and increases (in our case) 
quadratically with the particle momentum. In the limit 
77/s— 5-0 (i.e. for ideal fluid dynamics) 6f vanishes. 

Since VISH2+1 is a (24-l)-dimensional code with longi- 
tudinal boost invariance but UrQMD propagates particles 
in all three spatial dimensions, we extend the VISH2+1 
output analytically from midrapidity {y = 0) to non-zero 
momentum rapidities, using boost invariance. After hav- 
ing determined the transverse position and transverse 
momentum of a particle using a space-time rapidity in- 
tegrated version of Eq. (1) as weight, we sample its 
momentum rapidity y randomly within the finite range 
—3 < y < 3, with a sharp cutoff at its upper and lower 
ends. (We restricted the range to |?/| < 3 to minimize the 
excitation of strings in UrQMD.) Its space-time rapidity rjs 
(defining its longitudinal position) is then sampled ac- 
cording to the 7?s-dependence of Eq. (1) (see Eq. (Al) for 
an explicit expression). This results in an 77^ distribu- 
tion of the generated particles that is fiat near rjs^O 
but smeared at the edges around r]s = ±3 (see inset 
in Fig. 15 for illustration). UrQMD propagates the re- 
sulting particles in all three spatial directions but (due 
to the nearly boost-invariant input over the y- and r]s- 
range mentioned) preserves boost invariance of the fi- 
nal momentum distributions accurately within the region 

< 1.5. This allows comparison of the calculations with 
midrapidity data from the RHIC experiments without 
having to worry about edge effects from the rapidity cut- 
off in H2O. 

The default switching temperatures used in this article 
are T^w = 165 MeV for s95p-PCE and Tsw = 160 MeV for 
SM-EOS Q(CE). Here, 165 MeV is equal to the chemi- 
cal freeze-out temperature Tchem used in s95p-PCE, and 
160 MeV is the transition temperature from the mixed 
to hadronic phase in SM-EOS Q [10]. In Sec. IV we also 
use lower Tgw values in order to study the existence of a 
switching window for VISHNU. 

For each hydrodynamic run we use H2O to generate 
a large number of events® whose particles are then fur- 



Sufficient accuracy for the p^-integrated V2 is obtained with 
event samples ranging from 2,000 for 0-5% centrality to 72,000 
for 70-80% centrality. For Figs. 1-3 we used 90,000 events in 
order to obtain sufficient statistics for V2{pt) of protons out to 
PT = 2 GeV/c. 



ther propagated with UrQMD until all collisions cease and 
unstable resonances have decayed. In the purely hy- 
drodynamic simulations we run the converter on the fi- 
nal decoupling surface with temperature T^ec to gener- 
ate events of free-streaming hadrons with similar statis- 
tics. (We use rdoc = 100MeV with EOS s95p-PCE and 
Tdoc = 130 MeV with SM-EOS Q (CE), corresponding to 
similar decoupling energy densities [32].) In these runs 
we have no further hadronic collisions and only allow the 
unstable resonances to decay, and analyze the final state 
with the same tools as used for complete VISHNU runs. 



C. Microscopic hadronic transport (UrQMD) 

For the modeling of the hadronic phase, we use the 
Ultra-relativistic Quantum Molecular Dynamics (UrQMD) 
model [18], a microscopic hadronic transport model based 
on the Boltzmann equation. UrQMD, which was initially 
developed as an ab-initio model for the simulation of rela- 
tivistic heavy-ion collisions, is well suited for the descrip- 
tion of a hadron gas in and out of chemical equilibrium 
[30, 43-45] and has been successfully applied to previ- 
ous hybrid micro-|-macro transport approaches based on 
ideal fluid dynamics [20, 23, 46]. 

In UrQMD, the system evolves through a sequence of bi- 
nary collisions and 2 ^ N decays of mesons and baryons. 
The cross sections are assumed to be free vacuum cross 
sections and depend on the center of mass energy of 
the two colliding hadrons as well as on their flavor and 
quantum numbers. The UrQMD collision term contains 
49 different baryon species (including nucleon, delta and 
hyperon resonances with masses up to 2 GeV) and 25 
different meson species (including strange meson reso- 
nances), which are supplemented by their correspond- 
ing anti-particle and all isospin-projected states. Full 
baryon/antibaryon symmetry is included. For excita- 
tions with higher masses a string picture is used. All 
states listed can be produced in string decays, s-channel 
collisions or resonance decays. 

Tabulated or parameterized experimental cross sec- 
tions are used when available, resonance absorption and 
scattering is handled via the principle of detailed bal- 
ance. If no experimental information is available, the 
cross section is either calculated via an OBE model or 
via a modifled additive quark model which takes basic 
phase space properties into account. A detailed overview 
of the elementary cross sections and processes included 
in the UrQMD model is given elsewhere [18]. 

When modeling the hadronic phase subsequent to the 
decay of a thermalized QGP, we flnd that the relative 
momenta and cm. energies of the individual hadron- 
hadron interactions are rather small and therefore string 
excitations and decays arc strongly suppressed and occur 
rarely. The dominant forms of interactions encountered 
are elastic scattering and inelastic scattering through res- 
onance formation and decay. 

It is important to note that UrQMD makes no equilib- 
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rium assumptions and can therefore be utilized for the 
description of systems in and out of equihbrium. UrQMD 
will retain equilibrium if given an equilibrium initial con- 
dition with suitable boundary conditions. This char- 
acteristic allows for the construction of our hybrid ap- 
proach, since it ensures that UrQMD's response to the mi- 
croscopic hadronic configuration generated by H2O will 
initially mimic the response of the viscous hydrodynamic 
model. However, the real advantage of transitioning to a 
microscopic transport such as UrQMD lies in its ability to 
describe the evolution of systems out of equilibrium, e.g. 
prior to equilibration or during the break-up and freeze- 
out stage of the reaction when assumptions of chemical 
and/or kinetic equilibrium are no longer valid. 

By virtue of its microscopic nature, UrQMD takes the 
full local temperature and particle fugacity dependence 
of the hadronic viscosity into account [45], even though 
it may be challenging to quantify the exact value of the 
hadronic viscosity inherent in the UrQMD calculation of 
the hadronic phase (due to its dependence not only on 
temperature but also on multiple (non-equilibrium!) par- 
ticle fugacities). 



III. SPECTRA AND FLOW FROM 
HYDRODYNAMICS AND THE HYBRID MODEL 

A. EOS dependence: s95p-PCE vs. SM-EOS Q 

In this subsection we calculate transverse momentum 
spectra and elliptic flow for identified hadrons in the hy- 
brid approach and compare them to pure ideal and vis- 
cous hydrodynamic calculations, focusing in particular 
on the EOS and the Tdoc dependences of the calcula- 
tions. The main idea of this comparison is to determine 
whether the results of the hybrid model calculation can 
be reproduced by pure ideal or viscous hydrodynamics, 
or if specific features in the spectra and the transverse 
momentum or centrality dependence of V2 lead to clear 
discriminators between the two approaches. Unless oth- 
erwise noted, we use initial conditions and switching and 
decoupling parameters as described in Sec. II. With iden- 
tical initial conditions and identical transport coefficients 
r]/s in the hydrodynamic QGP stage, we ensure that the 
only difference between the VISH2+1 and VISHNU runs 
lies in the treatment of the hadronic stage. In the SM- 
EOS Q(CE) case this means that the pure hydrodynamic 
runs assume both thermal and chemical equilibrium in 
the hadronic phase, whereas the UrQMD component of 
VISHNU breaks both. When we use VISH2+1 with s95p- 
PCE, we ensure that the pure hydrodynamic and VISHNU 
hybrid simulations use the same non-equilibrium chem- 
ical composition in the hadronic stage; in this case the 
difference lies only in the assumption of (approximate) 
thermal equilibrium for the hydrodynamic runs whereas 
the UrQMD cascade in VISHNU allows the system to evolve 
far out of local thermal equilibrium, all the way to final 
decoupling. 



The pion and proton transverse momentum spectra 
shown in the left panels of Fig. 1 demonstrate good agree- 
ment between the pure hydrodynamic and hybrid VISHNU 
runs as long as EOS s95p-PCE is used, i.e. as long as we 
ensure that the hydrodynamic simulations correctly im- 
plement the non-equilibrium chemical evolution in the 
hadronic phase. With the correct PCE equation of state, 
the hydrodynamic pion spectra are almost insensitive to 
the choice of kinetic decoupling temperature Tdcc, agree- 
ing well with their counterparts from the hybrid model 
VISHNU in all cases. The hydrodynamic proton spec- 
tra become flatter as Tdoc decreases, due to build-up 
of additional radial flow, and the best agreement with 
the hybrid model is achieved for the lowest value shown 
in the graph (Tdoc = 100 MeV). This demonstrates that 
significant additional radial flow is generated during the 
hadronic rescattering stage [21, 47], but that for pions 
(whose momentum distributions react less strongly to ra- 
dial flow than the heavier protons) the radial flow and 
cooling effects on the spectral slope balance each other 
as we lower Tdec [35] . Figure 1 shows that all of the above 
statements hold irrespective of whether we assume zero 
or non-zero viscosity during the QGP stage: in both cases 
the transverse momentum spectra from the hybrid code 
VISHNU can be well represented by purely hydrodynamic 
simulations with Tdec = 100 MeV and rj/s values that do 
not change between the QGP and hadron gas stages. 

The right panels in Fig. 1 show that the same is 
not true when we use SM-EOS Q in the hydrodynamic 
code which assumes chemical equilibrium in the hadronic 
stage, contrary to the underlying microscopic UrQMD dy- 
namics. In this case agreement between the pion spectra 
from VISHNU and the purely hydrodynamic simulations 
requires immediate decoupling at Tgw = 160 MeV, but the 
corresponding purely hydrodynamic proton spectra are 
too steep because they lack the boost from the additional 
radial flow developed by UrQMD during the hadronic stage. 
Lowering Tdoc hi the pure hydrodynamic runs helps with 
the shape of the proton spectra, but quickly eats into 
the total proton yield (i.e. the normalization of the pro- 
ton spectra), due to annihilation with antiprotons, and 
simultaneously the pion spectra become too fiat when 
compared with VISHNU. 

We conclude that, at the level of single-particle pr 
spectra, it is possible to simulate the full microscopic 
dynamics of the UrQMD hadron cascade by viscous (or 
even ideal) fluid dynamics as long as a PGE EOS 
is used that correctly describes the non-equilibrium 
chemical composition in UrQMD, and one allows for the 
buildup of additional hadronic radial flow by setting a 
low dccouphng temperature Tdoc 100 MeV. With SM- 
EOS Q(CE), purely hydrodynamic simulations are un- 
able to reproduce the VISHNU spectra for any choice of 
Tdcc; for the convenience of the following academic com- 
parison of elliptic flow V2 within these two approaches, 
we choose for SM-EOS Q the "historical standard value" 
Tdcc = 130 MeV [48]. 

Research over the past few years has established that 
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FIG. 1. (Color online) Pion and proton pr-spectra for 6 = 7fm Au+Au, calculated from VISH2+1 and VISHNUwith s95p-PCE 
(left panels) and SM-EOS Q(CE) (right panels). Upper panels compare pure ideal fluid dynamics and VISHNU with ideal fluid 
input; lower panels compare pure viscous hydrodynamics and VISHNU with viscous hydrodynamic input, using t]/s = 0.08 in 
both cases during the fluid dynamic stage. Identical Glauber model initial conditions are used in all runs (eo = 30GeV/fm^ at 
To = 0.6fm/c). 



elliptic flow V2 is influenced by both hadronic dissipativc 
effects [22] and the chemical composition of the hadronic 
matter [34]. The reaction dynamics modeled by UrQMD 
contains both types of non-equilibrium effects. To isolate 
the effect of kinetic non-equilibrium we can compare re- 
sults from VISHNU with pure hydrodynamic calculations 
that model the same hadronic chemical non-equilibrium 
composition as UrQMD, by using EOS s95p-PCE. This is 
shown in Fig. 2(a), for both an ideal (77/5 = 0) and min- 
imally viscous 77/5 = 0.08 QGP fluid. For (77/s)qgp = 
one finds that in the hybrid model calculations V2 is 
suppressed by ^ 15% at pr — 2 GcV, as a result of the 
hadronic viscosity inherent in the UrQMD dynamics. For 
77/5 = 0.08 the difference between the purely hydrody- 
namic and hybrid model results is a bit less (~ 10%) since 
now the non-zero hydrodynamic viscosity already par- 
tially accounts for the viscous V2 suppression in UrQMD. 
Comparing the ideal and viscous pure hydrodynamic re- 
sults with each other we find ^ 20% viscous V2 suppres- 
sion a.t pt — 2 GeV for r]/s = 0.08, consistent with earlier 
results [8-10]. Fig. 2(a) demonstrates that evolution with 
the hybrid model suppresses V2 more strongly than vis- 
cous hydrodynamics alone with the same 77/s; a realistic 



microscopic description of dissipative hadron dynamics 
within a hybrid approach is therefore essential.^ 

Figure 2(b) shows that an incomplete separation of 
chemical and kinetic non-equilibrium effects can lead to 
misleading conclusions. Using a chemical equilibrium 
EOS (SM-EOS Q) in the hadronic phase and compar- 
ing a pure ideal fluid calculation with a VISHNU simu- 
lation with ideal fluid input (which is out of chemical 
equilibrium during most of the hadronic stage) happens 
to yield (at least for 6 = 7 fm in Au-I- Au^° ) almost iden- 
tical results for V2{pt)- Viscous suppression of V2 by 
kinetic non-equilibrium [22] almost exactly balances [21] 



^ We will see below that merely adjusting the specific shear vis- 
cosity r)/ s to larger values in the hadronic phase is not sufficient. 

" When compared with experiment, the hybrid model calculations 
exhibit a significantly improved centrality dependence compared 
to pure ideal fiuid dynamics. The almost perfect cancellation of 
kinetic and chemical non-equilibrium effects on 1^2 (pt) seen in 
Fig. 2(a) may therefore not work as well at other impact param- 
eters, although the tendency of the two effects to work against 
each other is generic [25]. 
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FIG. 2. (Color online) v-2.(j)t) for all hadrons from & = 7fm Au+Au collisions, calculated from VISH2+1 (lines without symbols) 
and VISHNU (lines with symbols). Dashed and solid lines are for ideal (77/5 = 0) and viscous (77/5 = 0.08) fluids during the 
hydrodynamic stage, with EOS s95p-PCE (a) and SM-EOSQ Q(CE) (b), respectively. Same initial conditions as in Fig. 1. 



the previously observed enhancement of V2(j)t) caused 
by chemical non-equilibrium in the hadronic phase [32- 
34]. The solid lines comparing pure viscous hydrodyna- 
mics and VISHNU with viscous fluid input demonstrate 
that this cancellation is accidental and no longer occurs 
when describing the QGP as a viscous fluid. In this case 
VISHNU gives much lower 1)2 (pt) than the pure viscous hy- 
drodynamic calculation, mainly caused by large negative 
contributions from the viscous correction 8j to the local 
distribution function on the switching surface which the 
subsequent UrQMD dynamics is unable to erase. Viscous 
hydrodynamics with constant 77/s in the hadronic phase, 
on the other hand, evolves back towards local thermal 
equilibrium such that on the final decoupling surface at 
Tdoc only small 5/ contributions remain. 

B. Hadronic contribution to elliptic flow 

Having established that identical chemical composi- 
tions in the pure hydrodynamic and hybrid model evolu- 
tions are essential for a meaningful comparison that aims 
to assess dissipative effects in the hadronic stage, we will 
from now on use the chemically frozen EOS s95p-PCE in 
all simulations. 

In Fig. 3 we compare the differential pion and pro- 
ton elliptic flow, V2{pt), from VISHNU (solid lines 
with symbols) with pure fluid dynamical calcula- 
tions, implementing kinetic freeze-out either directly at 
Tsv, = Tchcm = 165 MeV (dashed lines with symbols) or at 
Tdec = 100 MeV (solid lines without symbols). The curves 
in panel (a) and (b) assume ij/s — O and 7y/s = 0.08; re- 
spectively, during the fluid dynamic stage. We focus our 



This is, of course, an academic comparison since viscous hy- 
drodynamics with constant r]/s cannot consistently account for 
hadronic freeze-out. 



attention on the mass splitting between pions and pro- 
tons. The smallest mass splitting is found for immediate 
decoupling at Tgw (dashed lines). Hadronic rescattering 
in VISHNU (solid lines with symbols) increases the mass 
splitting at low pt, by pushing V2{pt) up for pions and 
towards larger pt for protons [47] . Both of these effects 
strengthen if we replace the microscopic hadronic rescat- 
tering cascade by macroscopic hydrodynamic evolution 
[47] . The depletion of proton V2 at small px is a conse- 
quence of additional radial flow buildup in the hadronic 
stage which pushes the heavy protons more efhciently 
than the light pions to larger tranverse momenta (see 
Fig. 1); the larger V2{pt) for low-p^ pfons reflects (at 
least partially) a simultaneous increase of the total mo- 
mentum anisotropy, in response to the remaining spatial 
fireball eccentricity that survives into the hadronic stage. 
These effects are qualitatively similar for ideal (Fig. 3a) 
and viscous fluids (Fig. 3b), although the larger 5f cor- 
rections at Tsw in viscous hydrodynamics lead to an ad- 
ditional downward shift of V2{pt) for both pions and pro- 
tons at large pt when hadronic rescattering is turned off 
(dashed lines). 



C. Viscous V2 suppression 

Figure 4 shows the net effect of the hadronic 
medium modifications of the transverse momentum spec- 
tra and differential elliptic flow on the pT-integrated 
charged hadron V2, as a function of collision centrality 
(parametrized by the number of participants A^part)- The 
smallest amount of V2 is obtained without any hadronic 
rescattering at all (dashed lines) , with an additional sup- 
pression of about 20% for the minimally viscous fluid 
(?7/s = 0.08) relative to the ideal fluid. The UrQMD mod- 
ule in the hybrid code VISHNU creates about 15% addi- 
tional V2 via hadronic rescattering, but a hydrodynamic 
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FIG. 3. (Color online) Differential elliptic flow D2(pt) for pions (light blue/magenta) and protons (black/purple) from fe = 7fm 
Au+Au collisions. Results from the hybrid model VISHNU(solid lines with solid symbols) are compared with purely hydrodynamic 
simulations for 77/3 = (a) and r]/s = 0.08 (b) for decoupling temperatures Tdec = 165 MeV (dashed with open symbols) and 
Tdec = 100 MeV (solid without symbols). All calculations use EOS s95p-PCE in the hydrodynamic stage. Initial conditions are 
the same as in Fig. 1. Please note the suppressed zeroes on the horizontal axes. 



description of the hadron gas, treating it either as an 
ideal or a minimally viscous fluid, generates a much larger 
hadronic contribution to V2 ■ This reflects the fact that in 
the hadronic stage the fireball is still out-of-plane elon- 
gated, and demonstrates that even a viscous fluid with 
?7/s = 0.08, but especially an ideal fluid is much more 
efficient than UrQMD in converting this residual fireball 
eccentricity into additional elliptic flow. 

To quantify the suppression of integrated elliptic flow 
by viscous QGP and hadronic dissipation effects we de- 
fine the ratio 



where A and B denote two different dynamical evolu- 
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FIG. 4. (Color online) Gentrality dependence of the py- 
integrated elliptic flow of all hadrons, for the same parameters 
as in Fig. 3. 



tion models. Figure 5a shows that viscous suppression 
effects are generically larger in peripheral than in cen- 
tral collisions, due to the smaller fireball sizes. The solid 
line in panel (a) confirms the naive expectation that the 
strongest suppression (here 30-50%, depending on colli- 
sion centrality) should be seen when comparing VISHNU 
with viscous fluid input to a purely hydrodynamic, ideal 
fluid simulation. Among the combinations studied in 
Fig. 4, this case maximizes the suppression of elliptic flow 
by combining (in the VISHNU simulation) viscous effects 
in the QGP with large dissipative effects in the hadronic 
stage. Replacing the UrQMD part of the evolution by vis- 
cous hydrodynamics with 77/5 = 0.08 (dashed line) yields 
only about half of the viscous vi suppression observed in 
VISHNU. The smallest suppression ratio (~ 10 - 15%) is 
found between VISHNU simulations with ideal and viscous 
fiuid input at Tg^ (dot-dashed line in Fig. 5a). 

These observations indicate that, in Au+Au collisions 
at RHIC with (?7/s)qgp = , hadronic dissipative 
effects play a larger role for the finally observed v-i than 
QGP viscosity. This conclusion is reinforced by Fig. 5b 
with shows (depending on centrality) 20-40% V2 suppres- 
sion just from hadronic dissipation. The VISH2+1 and 
VISHNU comparison runs for 77/s = l/(47r) show relatively 
less suppression of vi by hadronic dissipation in UrQMD 
since the non-zero 77/s already suppresses hadronic vi 
generation in the purely hydrodynamic run. 

Obviously, the hadronic rescattering stage described 
by UrQMD is much more dissipative than both an ideal 
(?;/s = 0) and a "perfect" (i.e. minimally viscous, ^ = ^) 
fluid. This raises the question whether we could perhaps 
simulate the hadronic rescattering cascade hydrodynam- 
ically by making the fluid more viscous on the hadronic 
side of the quark-hadron phase tansition. The answer to 
this question is explored in the following section. 
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IV. VISHNU WITH LOWER SWITCHING 
TEMPERATURES AND A TEMPERATURE 
DEPENDENT HADRONIC {'n/s){T) 

Hydro+cascade hybrid approaches that use in their hy- 
drodynamic modules an equation of state that assumes 
chemical equilibrium in the hadron resonance gas have 
little choice where to switch from the macroscopic fluid 
dynamic to the microscopic Boltzmann picture: If they 
want to correctly reproduce the experimentally measured 
final hadron yields which reflect chemical frcczc-out at 
Tchem — 165 McV, the switching must be done at that 
temperature. With the chemically frozen EOS s95p-PCE 
we can also select lower switching temperatures and, in 
doing so, explore the existence of a "switching window" 
within which, using appropriately adjusted transport co- 
efficients for the hydrodynamic evolution, both macro- 
scopic and microscopic descriptions can be used inter- 
changeably, without affecting the final outcome. 

To judge the validity of using hydrodynamics to emu- 
late microscopic UrQMD dynamics we check the final pion 
and proton transverse momentum spectra and their dif- 
ferential elliptic flow V2 {pt ) ■ Wc do so for a fixed impact 
parameter 6 = 7 fm which for the Au+Au collision sys- 
tem is known to also provide a fair representation of the 
spectra and elliptic flow from minimum bias collisions. 

The most sensitive and robust [12] observable that will 
dictate our (temperature dependent) choice of rj/s in the 
hadronic phase is the integrated elliptic flow V2- Strictly 
speaking, it is controlled by both the shear viscosity and 
the associated microscopic relaxation time r,r that con- 
trols the speed with which the shear pressure tensor tt'^" 
approaches its Navier-Stokes limit [10]. In kinetic trans- 
port theory, both 77 and involve integrals over the 
same collision kernel and are typically proportional to 
each other [49]. In the QGP, where r]/s is small, is 
therefore short of order 0.2 fm, leading to rapid memory 
loss and insensitivity of the developing elliptic flow to ini- 
tial conditions for the viscous pressure [10]. If, however, 
r]/s becomes large in the hadronic phase [45, 50-52], one 



should expect the corresponding relaxation time to grow 
similarly, with a constant of proportionality between the 
dimensionless combinations 77/s and Ttt^ that may be 
different and possibly larger in the hadronic than in the 
QGP phase. 

Unfortunately, with our present code setup an inde- 
pendent determination of the temperature dependence 
of rj/s and Tt^^ in the hadronic stage (from a comparison 
of VISHNU simulations with different Tgw values) is too 
time consuming and thus not practical. We have there- 
fore focused our attention on the extraction of (r]/s)(T), 
holding the relation Ttt^ = 2>ri/ s fixed at all temperatures. 
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A. Effective {t]/s){T) from UrQMD with VISH2+1 input 

The green dashed line in Fig. 6 shows that, if we cou- 
ple VISH2+1 with constant specific entropy ?//s = 0.08 
to UrQMD at lower and lower switching temperatures, 
we obtain larger and larger values for the total pion 
1)2. Obviously, the hydrodynamic evolution between 
Tchcm = 165 McV and Tsw < Tchem generates more ad- 
ditional elliptic flow than propagation of the hadrons 
with UrQMD during the same temperature interval. 
This suggests that UrQMD in the temperature region 
^sw <T < Tchom has a larger effective shear viscosity 
r]/s than 0.08. By lowering the switching temperature 
from the starting value T^w = 165 MeV in small steps 
AT and adjusting, step by step, rj/s in the interval 
[165MeV-nAr, 165 MeV-(n-l)AT] such that no ad- 
ditional pion V2 is generated when we replace UrQMD evo- 
lution in this temperature interval by VISH2+1 with this 
adjusted rj/s value, we arrive at the effective tempera- 
ture dependent (r^/s)^^-* (T) shown in the inset of Fig. 6) 
that ensures a constant pion elliptic flow V2 that is in- 
dependent of the switching temperature (flat horizontal 
line in Fig. 6). Note that this is a time consuming and 
labor-intensive iterative procedure that cannot be short- 
circuited since the effective 77/s extracted for the n"^ in- 
terval turns out to depend on the previously determined 
effective {ri/s)^^\T) values for smaller n values, i.e. at 
higher temperatures. 

B. Testing the effective [r] / s)^^'> {T) from UrQMD 

Starting at the QGP input value of 7//s = 0.08, the 
thus extracted effective {r]/s)'^^^T) increases by almost a 
factor 5 between T = Tchom = 165 MeV and T = 100 MeV, 
growing roughly linearly with decreasing temperature. 



Figure 7 shows that this has very little effect on the 
pion and proton py-distributions: Whether one uses con- 
stant 77/5 = 0.08 (panel (a)) or the temperature depen- 
dent {rj/sY^^T) from Fig. 6 (panel (b)), the transverse 
momentum spectra for both pions and protons exhibit 
very little dependence on the switching temperature Tgw 
Just as for the pure hydrodynamic calculations shown in 
Fig. 1 , lower switching temperatures lead to slightly flat- 
ter spectra for protons (where the hadronic buildup of 
additional radial flow overwhelms the cooling effect) and 
to slightly steeper spectra for pions (for which the cooling 
effect dominates) . Very little of the difference in the final 
spectra that can be attributed to the different hadronic 
viscosities used in the two panels of Fig. 4. 



«_« vise, hydro + UrQMD; T =165MeV 
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FIG. 8. (Color online) Differential V2{pt) for pions and pro- 
tons in b = 7 fm Au-|- An collisions, calculated by matching vis- 
cous hydrodynamics with temperature dependent (r7/s)'^'(r) 
(see Fig. 6) to UrQMD at different switching temperatures Tsw . 

For the elliptic flow the temperature dependence mat- 
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ters, as Figure 6 clearly demonstrates. However, Fig- 
ure 8 shows that, when using the "correct" temperature 
dependent {?] / s)'^^\T) in the VISH2+1 module, VISHNU 
produces not only Tgw-independent integrated pion ellip- 
tic flow (by construction), but also approximately Tgw- 
independent pT- differential elliptic flow V2{pt) for both 
pions and protons. Still, when looking at the details 
one observes a difference between the solid and dashed 
brown lines for protons which shows that decoupling 
VISH2+1 at Tdcc — 100 MeV into non-interacting particles 
is not the same as switching from VISH2+1 to UrQMD at 
Tsw = 100MeV and letting UrQMD do the kinetic freeze- 
out: there obviously is some rcscattcring of protons in 
UrQMD occurring at temperatures below 100 MeV, caused 
by "pion wind" [53] , that moves the protons to larger pT 
and depletes the proton elliptic flow at low px- Since the 
solid brown line, corresponding to hydrodynamic evolu- 
tion down to 100 MeV followed by UrQMD freeze-out, is dif- 
ferent from the solid red line with triangles, which corre- 
sponds to a switch from hydrodynamics to UrQMD already 
at 165 MeV, we see that even with the "correct" tempera- 
ture dependent (Ty/s)*-^-* (T) hydrodynamic evolution with 
VISH2+1 differs somewhat from UrQMD. However, the dif- 
ference is small, and for pions its consequences are almost 
negligible. 

Figures 9 and 10 show that the same {rj / s)^^") (T) ex- 
tracted in Fig. 6 from simulations with optical Glauber 
model initial conditions also works for CGC motivated 
initial conditions that are obtained by averaging over 
many fluctuating initial entropy density profiles (see 
Sec. HA) and thus account for the large effects from 
event-by-event fluctuations on the average initial source 
eccentricity in very central and very peripheral events. 
For all three centrality classes shown in these two figures 
we see that using the temperature dependent (77/5)'-^'' (T) 
extracted from Fig. 6 for the hydrodynamic evolution 
yields values for the integrated pion elliptic flow V2 
(Fig. 9) and shapes for the pt-diffcrcntial pion elliptic 
flow V2{pt) (Fig- 10) that arc independent of the switch- 
ing temperature Tgw down to Tsw = 120MeV (we did 
not probe any lower). This is not trivial since the ini- 
tial source eccentricities (and thus the produced elliptic 
flows) for the centralities shown here are 20-50% larger 
than the optical Glauber ones. 

In the most peripheral bin (60-70%) things begin to 
break down around Tgw = 120 MeV. The increase of inte- 
grated V2 and of V2{pt) at low px seen in Figs. 9 and 10 
for Tsvi, = 120MeV indicates stronger dissipativc effects 
in UrQMD at low T than captured by effective (?7/s)(^'(T) 
from Fig. 6. In peripheral collisions (i.e. for small 
flrcballs) VISH2+1 apparently can no longer accurately 
mimic the microscopic kinetic evolution of UrQMD at low 
temperatures Tsw ~ 120 MeV and should therefore not be 
used under such conditions. 

We can summarize the findings of this subsection 
in the following statement: By using the temperature- 
dependent [ri/s)^'^\T) extracted in Sec. IV A from the 
integrated pion elliptic flow by comparing microscopic 
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and macroscopic simulations of the hadron gas phase 
below Tchem, thc macroscopic evolution code VISH2+1 
provides a fair description of the microscopic UrQMD dy- 
namics down to temperatures around 120 MeV. When 
used in hybrid inode together with UrQMD (as imple- 
mented in VISHNU), it yields transverse inomentum dis- 
tributions of pions and protons (including their ellip- 
tic flow V2(pt)) that are, to good approximation, in- 
dependent of the switching temperature in the win- 
dow 120 MeV < Tsw < 165 MeV. This raises, however, two 
questions: 

• Is (jy/s)*-^-* (T) a genuine medium property of the 
hadronic matter described by UrQMD? If it is, it 
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should be independent of the viscous hydrodynamic 
input into UrQMD. So far we have only shown results 
where the UrQMD input was obtained from hydro- 
dynamic simulations of the earlier QGP evolution 
with a single viscosity value (?7/s)qgp = 0.08. If the 
QGP has a larger or smaller viscosity, will UrQMD 
continue to propagate the correspondingly modified 
hydrodynamic input at Tgw as if it were a viscous 
fluid with temperature dependent {ij / s)^^'' {T)7 

• The {r] / s)'^^^ (T) values shown in the inset of Fig. 6 
are much smaller (especially in the region just be- 
low Tchem ~ Tc) than those computed in Ref. [45] 
using the Kubo formula. The Kubo formula eval- 
uates a thermal equilibrium ensemble expectation 
value of a spectral density operator, whereas the 
procedure of Sec. IV A extracts the transport co- 
efficient from a comparison of the dynamical out- 
put produced by two different evolution models in 
a rapidly expanding expanding medium. What is 
the reason for these different shear viscosities? 

In the following subsection we describe a study that 
yields at least partial answers to these questions. 



C. Non-universality of the effective UrQMD viscosity 

To address the first question we redid the analy- 
sis of Sec. IV A using a twice larger QGP viscosity, 
(^/•s)qgp = 0.16. This means that the UrQMD module of 
VISHNU is initialized with somewhat different density and 
flow profiles and, in particular, a shear viscous pressure 
tensor tt^'' ps 2r]a^'^ that is roughly twice as large as be- 
fore when we used (?7/s)qgp = 0.08.^^ Going through the 
same iterative procedure as before of adjusting the spe- 
cific shear viscosity r]/s temperature interval by temper- 
ature interval such that it generates in VISH2+1 exactly 
the same amount of total pion V2 as UrQMD, we obtain 
the upper set of solid green dots in Fig. 11, labeled as 
{rj / sY^^T) . It is obviously different from and larger than 
[r] / sY^'> (T) . In particular, at the highest switching tem- 
perature Tsw = 165MeV, (77/s)(^'(T) starts out close to 
the value of 0.16 used in the evolution of the QGP stage, 
whereas {r]/s)^^^{T) starts out almost half as small, with 
a magnitude close to the QGP value of 0.08 used in the 
first extraction. 

The "effective UrQMD viscosity" {'q/s){T) extracted by 
this procedure apparently "remembers" the prior QGP 
history of the fireball and its transport properties. This 



Here cr'"' = V^^'m''^ is the velocity shear tensor, see [10] for de- 
tails. Since rj/s is small in the QGP, the associated relaxation 
time Ttt is much shorter than the inverse of the local expansion 
rate on the switching surface, hence vr'"' does not stray far from 
its Navier-Stokes value tt^"^ = 2»;cr'"^ on that surface. 



means that it does not describe an intrinsic medium prop- 
erty of the hadron resonance gas in UrQMD. Further in- 
sight into this puzzle is gained by looking at the pur- 
ple dashed line and open squares in Fig. 11, labeled by 
{rj / sY^\T): Here we initialize UrQMD with viscous hydro- 
dynamic output from a run with (77/s)qgp = 0.16, but 
assume that just before hadronization, at a temperature 
of 180 MeV, the QGP viscosity suddenly drops to half 
that value, (77/s)qgp = 0.08. In this case, the extracted 
"effective UrQMD viscosity" (yy/s)'^'^^ (T) comes out identi- 
cal to (7?/s)(i)(T), i.e. as if the QGP phase had evolved 
with minimal shear viscosity throughout its life, and not 
just at the very end of its history after its temperature 
dropped below 180 MeV. 

Evolving the QGP medium from identical initial con- 
ditions with a twice larger t]/s value of 0.16 leads to 
noticeably different density and flow profiles on the 
Tsw = 165 MeV switching surface compared to those for 
77/s = 0.08. These differences cannot be undone by chang- 
ing rj/s back to 0.08 within 15 MeV of the critical temper- 
ature. The equality of (r;/s)(^'(r) and {ri/sY^\T) shows 
that they are not controlled by the energy density and 
flow velocity profiles on the switching surface (which are 
different for the solid red and dashed purple lines), but 
instead by the shear pressure tensor. Due to the short 
relaxation times r^r in the QGP phase, we can approxi- 
mate TT^'' « 2770-^'^, so a factor of 2 difference in 77/s leads 
to roughly a factor 2 difference in tt^'' on the switch- 
ing surface. (The differences in the entropy density and 
velocity shear tensor, although noticeable, are not any- 
where close to a factor 2.) At this level of precision, 
the hydrodynamic output corresponding to the green 
solid line (from a QGP evolving with (77/s)qgp = 0.16 
all the way to Tgw = 165 MeV) features a viscous pres- 
sure tt''" that is twice as big as that for the red solid line 
((77/s)qgp = 0.08), whereas the one corresponding to the 
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FIG. 11. (Color online) Temperature dependent effective 
{rj/s){T) extracted as in Fig. 6 from viscous hydrodynamic 
emulations of UrQMD dynamics, for three different viscous fluid 
initializations of the UrQMD stage (see text for discussion). 
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purple dashed line is roughly equal to that of the red 
solid line. The "effective UrQMD viscosities" {rj / s)'''^\T) 
seem to reflect and remember these relationships. 

We can understand this by recalling the kinetic theory 

relation tt^'" = / ^p^'' p "'^ Sf between the viscous pres- 
sure tensor and the deviation Sf of the phase-space dis- 
tribution function from local thermal equilibrium. The 
value of TT^'^ on the switching surface controls the magni- 
tude of the non-equilibrium contribution 6f ~ p^p'^T^^iv 
to the hadron momentum distributions sampled by the 
hydro-to-micro converter H2O which generates the UrQMD 
input. Apparently, collisions and other interactions in 
UrQMD arc too infrequent or inefficient in transferring 
momenta to quickly relax these non-equilibrium distri- 
butions. As long as the deviations 5f persist in a form 
close to their initial values, they contribute a term tt'"^ 
to the UrQMD energy-momentum tensor that, if one writes 
it in hydrodynamic language as tt^" = 2ria^'^ , requires a 
value of r]/s that is close to the value used in the gen- 
eration of the UrQMD input with H2O. As UrQMD evolves 
the momentum distributions, tt^"^ evolves, too, and when 
we fit the VISH2+1 transport properties to those of the 
UrQMD cascade assuming short relaxation times, we end 
up fitting this evolving tt^'^ instead of extracting a gen- 
uine transport coefficient for UrQMD. 

The (r;/s)(")(T) curves shown in Fig. 11 exhibit a ten- 
dency to approach each other at lower temperatures. 
This may suggest a loss of memory of the initial vis- 
cous pressure components on a time scale comparable 
to the cooling time needed to cool the system from 
rchem = 165MeV to somewhere around 100 MeV. This 
time is longer than the one assumed in our viscous hy- 
drodynamic simulations of the UrQMD stage, Tr^r = 377/5. 
Indeed, there are other indications [54] that UrQMD may 
have a much longer relaxation time than given by the 
classical kinetic theory result for massless Boltzmann 
particles, TTr = 6ri/{sT) [55]. For the time being, further 
exploration of this issue will have to wait. 



V. DISCUSSION AND CONCLUSIONS 

In this article we presented VISHNU, a new hybrid ap- 
proach for the bulk evolution of viscous QCD matter 
created in ultrarelativistic heavy-ion collisions that com- 
bines a macroscopic viscous hydrodynamic description of 
the early (dense) QGP phase with a microscopic Boltz- 
mann cascade for the late (dilute) hadron resonance gas 
stage. The model merges the economy of a macroscopic 
description for the first 4— lOfm/c (depending on cen- 
trality) after thermalization, when the matter is close 
to local thermal equilibrium, with the precision of a mi- 
croscopic approach for the late hadronic stage when, af- 
ter hadronization of the QGP, the mean free paths for 
the hadronic medium constituents increase rapidly, the 
matter moves farther and farther away from local equi- 
librium, and finally decouples into non-interacting, free- 
streaming particles. 



The hydrodynamic and Boltzmann cascade compo- 
nents of VISHNU are connected via the hydro-to-micro 
converter H2O, a Monte Carlo event generator that sam- 
ples the Cooper-Frye phase-space distributions from the 
viscous hydrodynamic output along a switching surface 
of constant temperature T^w and injects the resulting par- 
ticles into UrQMD for further microscopic propagation. 

By using a microscopic cascade approach for the late 
hadronic stage in heavy-ion collisions, the complex mi- 
croscopic dynamics of chemical and kinetic freeze-out 
is accounted for without extraneous parameters, remov- 
ing a major source of uncertainty in connecting final 
hadron spectra with transport properties of the initial 
QGP phase. With the help of VISHNU a rehable extrac- 
tion of the QGP shear viscosity (77/s)qgp from experi- 
mental RHIC data, with good theoretical control of the 
associated uncertainties, has now become possible [12]. 

In default mode, VISHNU is used with switching 
temperature Tgw = T'chem = 165 McV. Since the experi- 
mentally determined chemical freeze-out temperature 
Tchem ~ 165 MeV [28] approximately agrees with the 
best presently available theoretical estimate for the 
(pseudo)critical temperature for the chiral phase tran- 
sition in QCD [26, 27], this is the highest temperature 
for which a microscopic description in terms of colliding 
hadrons with vacuum masses and decay widths makes 
sense, i.e. the highest temperature for which UrQMD can 
be reliably used. In default mode the system is therefore 
described by a chemical equilibrium EOS in the hydro- 
dynamic stage, followed by a microscopic stage in which 
first chemical and then kinetic equilibrium are broken au- 
tomatically by the microscopic scattering dynamics. In 
this mode no chemically frozen EOS is needed. 

The main thrust of the present work has been to in- 
vestigate whether VISHNU can be used with lower switch- 
ing temperatures. Since the microscopic evolution of 
the hadron resonance gas via UrQMD is numerically much 
more costly than a macroscopic hydrodynamic descrip- 
tion, one would like to use hydrodynamics as long as 
possible and switch to UrQMD only when the macroscopic 
approach is no longer reliable. We showed that in that 
case using a chemically frozen hadron resonance gas EOS 
in the hydrodynamic evolution below Tchcm is compul- 
sory. Without such a realistic EOS, which takes into ac- 
count that at temperatures below Tchem the stable par- 
ticle yields must not change anymore, one cannot find 
a single value for the kinetic decoupling temperature 
in VISH2+1 that simultaneously reproduces the micro- 
scopically evolved pion and proton transverse momen- 
tum spectra from VISHNU. For the differential elliptic flow 
V2{pt) using a chemical equilibrium hadronic EOS can 
lead to quite misleading results. 

Using the chemically frozen EOS s95p-PCE [24, 25] one 
finds that, with the judicious choice Tdcc — 100 MeV for 
the kinetic decoupling temperature, the pion and proton 
spectra from VISHNU can be reasonably well reproduced 
by a purely hydrodynamic calculation with VISH2+1, 
without even a need for changing the specific shear vis- 
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cosity T]/ s between the QGP and HRG stages. The same 
does not hold, however, for the elhptic flow: since at 
RHIC energies the fireball is still out-of-plane elongated 
when its matter enters the HRG phase, the momen- 
tum anisotropy continues to grow during the hadronic 
stage, and UrQMD is less efficient in converting the resid- 
ual source eccentricity into elliptic flow than VISH2+1 if 
one uses the same minimal t//s in both phases. Since 
we also showed that, under RHIC conditions, the total 
suppression of elliptic flow below its ideal fluid limit is 
dominated by the hadronic stage (Fig. 5), this demon- 
strates unequivocally that, if one wants to avoid or mini- 
mize the cost of using a microscopic approach like UrQMD, 
great care must be taken to use the correct transport 
properties of the hadronic phase in VISH2+1 below Tc. 

We therefore tried to improve the hydrodynamic de- 
scription of the HRG by increasing 77/s in the hadronic 
phase. To extract the temperature dependence of 77/s 
in the hadronic phase as described by UrQMD, we low- 
ered Tsw in small intervals and adjusted 77/s in VISH2+1 
in each newly covered temperature interval such that 
VISH2+1 produced exactly as much additional V2 for pi- 
ous as VISHNU resp. UrQMD did. We found a function that 
we called (?7/s)(^)(T) that started at Tgw^Tchcm close 
to the value used in the QGP phase and then increased 
dramatically, by almost a factor 5, as we lowered Tg^i to 
100 MeV. Surprisingly, this {r^/s)^^\T) was much smaller 
than the corresponding values previously extracted from 
UrQMD using the Kubo formula [45]. 

With this (77/s)(i)(T), VISH2+1 is able to reproduce 
quite well the transverse momentum spectra as well as 
not only the integrated, but also the pT-differential el- 
liptic flow of pious and protons calculated with the full 
hybrid code VISHNU (Figs. 7 and 8). It continues to do so 
if we replace the optical Glauber model initial conditions 
in VISH2+1 by an ensemble average of fluctuating initial 
conditions from the fKLN model (Figs. 9 and 10). But 
it fails badly when we replace the hydrodynamic input 
into UrQMD by one that was calculated with a different 
(77/s)qgp- In fact, the effective hadronic {rj/s){T) ex- 
tracted from UrQMD by demanding that both VISH2+1 
and VISHNU produce the same total pion elliptic flow ap- 
pears to track the shear viscosity of the preceding QGP 
phase: lower values of (77/s)qgp produce lower values of 
{•q/ s){T) below Tc, and vice versa. Clearly, the extracted 
hadronic {rj/ s){T) is not an intrinsic medium property of 
the hadron resonance gas in UrQMD, but a parameter that 
preserves some memory of the QGP transport properties 
(in fact, only of the transport properties at temperatures 
just above Tchom, see Fig. 11). 

Through a sequence of cross-checks we convinced our- 
selves that our method to extract an effective {rj/ s){T) for 
the hadronic matter in UrQMD failed because in UrQMD the 
viscous corrections in the energy-momentum tensor relax 
to their Navicr-Stokes values much more slowly than as- 
sumed in VISH2+1. As a consequence, the extracted 77/s 
fitted the value of the viscous pressure tensor tt'"^ inher- 
ited by UrQMD from VISH2+1 (and thus characterized by 



the 77/s value in the preceding QGP phase) instead of the 
shear viscosity characteristic of UrQMD matter. 

An independent extraction of both the relaxation time 
and the shear viscosity from UrQMD is presently beyond 
our technical means. Our results suggest that it may be 
possible to achieve a good emulation of UrQMD dynam- 
ics with viscous hydrodynamics, by using VISH2+1 with 
larger specific shear viscosities 77/s combined with larger 
relaxation times in the hadronic phase than those as- 
sumed in the present work (which were constrained by 
the relation Tr^r = 3 77/s). We are skeptic, however, about 
the prospects for finding a pair of functions {rj/ s){T) and 
Ttt,{T) that work universally, i.e. that yield good vis- 
cous hydrodynamic emulations of UrQMD independent of 
the initial input into UrQMD (in particular, independent 
of the initial values of the viscous terms in T^^'^ that 
are to be further evolved with UrQMD). The reason for 
our skepticism is that large relaxation times r^r indicate 
a fundamental breakdown of the viscous hydrodynamic 
framework: They appear as second-order corrections to 
ideal fluid dynamics in a systematic gradient expansion, 
and if they lead to large excursions of the viscous pres- 
sure TT^'^ away from its first-order Navier-Stokes value 
tt'^" = 277ct'"^, this indicates that the gradient expansion 
is no longer converging. 

With this state of present knowledge we conclude that 
there exists no "switching window" of temperatures be- 
low Tchom for safely switching from VISH2+1 to UrQMD, 
i.e. there is no temperature interval below Tchem iu which 
VISH2+1 and UrQMD can both be used equally well to de- 
scribe the evolution of the expanding medium created in 
relativistic heavy-ion collisions. Quantitative compari- 
sons with experimental data, with the goal of extract- 
ing from measured hadron spectra precise information of 
QGP transport properties, will therefore necessarily re- 
quire the use of a hybrid code like VISHNU in which the 
dynamics of the late hadronic stage (at all temperatures 
below Tchem) is described microscopically. 
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FIG. 12. (Color online) pr spectra for central Au+Au collisions and V2{pt) for non-central Au+Au collisions at b=7 fm, from 
VISH2+1 and H2O. 



Appendix A: Verification of the VISH2+1 to UrQMD 
converter 

In VISHNU the connection between VISH2+1 and UrQMD 
is realized by the Monte-Carlo particle generator H2O, 
which is briefly described in Sec. II B. Tliis appendix de- 
tails the verification of H2O against the analytical Cooper- 
Frye formula on which it is based. 

H2O randomly generates an ensemble of particles in 
momentum and position space based on the differential 
Cooper-Frye formula, Eq. (1), for each collision event, 
which can subsequently be used as initial configuration 
for UrQMD. Here, we directly perform a statistical analysis 
of sufficiently many such initial configurations obtained 
from H2O to generate smooth particle spectra and elliptic 
flow curves, which we then compare to a direct numer- 
ical integration of Eq. (1) along the freeze-out surface 
generated by VISH2+1. 

Figure 12 shows p^-spectra for pions, protons and 
kaons in central (a) and their differential V2{pt) hi non- 
central (b) Au+Au collisions, obtained from VISH2+1 di- 
rectly and via H2O, respectively. To specifically test the 
"viscous part" of the H2O event generator, we select a 
constant specific shear viscosity 77/s = 0.08. With a suffi- 
ciently large number of events, the statistical results from 
H2O exactly reproduce the hydrodynamic spectra and V2 
directly obtained from VISH2+1. 

The above hadron spectra are generated on a decou- 
pling hypersurface with Tdcc = 130 MeV.^'^ In Fig. 13, we 
show the 2-dimensional r—r freeze-out hypersurface for 
central Au+Au collisions. The black line is the hydro- 
dynamic freeze-out hypersurface obtained from VISH2+1, 
and the red symbols arc the pion emission points taken 
from a single H2O event, all of which fall on the VISH2+1 



freeze-out hypersurface. The small deviations which can 
be seen are due to the finite position resolution in H2O, 
which randomly samples particle positions on the freeze- 
out hypersurface within Ar = 0.1fm. Fig. 13 also shows 
that only few particles come from the very early stage 
of the fireball evolution, while most of the particles are 
emitted during the middle and late stages when the flow 
velocity is fully developed. This is illustrated in the inset 
of Fig. 13, which shows the emission rate for pions as 
a function of time (along the freeze-out hypersurface). 
Again the statistical results from H2O are in excellent 
agreement with the hydrodynamic results from VISH2+1 
showing that H2O correctly reproduces the particle emis- 
sion rates of VISH2+1 along the freeze-out hypersurface. 

Momenta and positions of the produced particles are 
correlated through Eq. (1). For example, the equilibrium 
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We here selected the "default" Tdoc = 130MeV for EOS-Q, but 
the quality of the code verification results docs not depend on 
the chosen decoupling/switching temperature. 



FIG. 13. (Color online) VISH2+1: hydrodynamic freeze-out 
surface for central (6 = 0) Au+Au collisions with decoupling 
temperature Tdcc ~ 130 MeV. H2O: pion emission points along 
the freeze-out surface from a single event. Inset: pion emission 
rate along the freeze-out surface. 
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— VISH2+1: hydrodynamic results 
* HjO: statistical results from 1000 events 
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FIG. 14. (Color online) Left: Pion emission probabilities as 
a function of emission angle A0= |0p — 0„| from VISH2+1 and 
H2O. 



distribution / for fi = can be explicitly written as 

1 



fcq ~ 



g7i[mT cosh(j/-77a)-pTfj_ cos(0p-</)„)]/T _|_ 2' 



(Al) 



This shows directly the correlation between the trans- 
verse flow velocity v± and transverse momentum px 
of the particles, as well as the correlation between 
monrentum rapidity y and space rapidity r]s- (Here 
(j)p = arctan^px / Py) is the angle of the particle transverse 
momentum and 0u = aTctan{vx/vy) is the angle of the 
transverse flow velocity.) 

In fact, the above form of the distribution function 
leads to enhanced particle production along the flow ve- 
locity direction, but suppresses particle production in the 
opposite direction. This is shown in Fig. 14 which shows 
the particle production rate as a function of the relative 
emission angle A(j)= \(f>p—(j)y\. Both the hydrodynamic 
result from VISH2+1 as well as the statistical result from 
H2O show that AN/ Ac/) reaches a peak when A(f> is zero 
(particle emission along the flow velocity direction), then 
rapidly decreases with increasing A(j), and reaches a min- 
imum for A(^=180°. The excellent agreement between 
the VISH2+1 and H2D results show that H2O correctly de- 
scribes the particle momentum and flow velocity correla- 
tions encoded in Eq. (1). 

The Bjorken approximation in (2-l-l)-dimensional vis- 
cous hydrodynamics leads to a uniform particle density 



as a function of momentum-space rapidity and likewise a 
uniform particle density as a function of space-time ra- 
pidity. Eq. (1) (together with Eq. (Al)) shows the corre- 
lation between momentum-space and space-time rapidity, 
which prevents H2O from generating particles with inde- 
pendent momentum-space and space-time rapidities. In 
practice, H2O first randomly generates the momentum- 
space rapidity y for each hadron within a finite pre- 
defined range (for example between -3 and -1-3) and then 
samples the space-time rapidity rjs through Eq. (1). The 
finite range of y, together with the y—r/s correlation, leads 
to a decrease of dN/ dijs in the boundary region of ris as 
shown in the inset of Fig. 15. 

Fig. 15 shows the correlation between y aird rjs for pi- 
ous. In excess of 90% of pions are produced at a space- 
time rapidity which lies within one unit of its momentum- 
space rapidity. In other words, most of the particles with 
forward (backward) momentum-space rapidity are gener- 
ated in forward (backward) space-time rapidity regions. 
Fig. 15 also shows that the hydrodynamic results from 
VISH2+1 and statistical results from H2O agree very well 
with each other, indicating that H2O correctly describes 
the momentum-space and space-time rapidity correlation 
in particle production. 
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— VISH2+1: hydrodynamic results 
* HjO: statistical results for 5000 events 



AY=lTi^-yl 
"T|^: space rapidity^ 
y: momentum rapidit>^ 



1 ' 1 ' 1 ' 1 ' 1 _ 










•— • dN/dy 






« dN/dti^ 




1 ^/s = 0.08 




- / =130 MeV \ - 

J dec V 

< 1 , 1 < 1 < 







FIG. 15. (Color online) Left: Pion emission probabilities as 
a function of emission rapidity AY= |?7s— J/| for VISH2+1 and 
H2O. Inset: Pion production as a function of momentum ra- 
pidity y and space-time rapidity rjs, from H2O. 
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